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ABSTRACT 

We report on Chandra ACIS imaging of the pulsar wind nebula (PWN) of the young Vela-like 
PSR B1706— 44, which shows the now common pattern of an equatorial wind and polar jets. 
The structure is particularly rich, showing a relativistically boosted termination shock, jets with 
strong confinement, a surrounding radio/X-ray PWN and evidence for a quasi-static 'bubble 
nebula'. The structures trace the pulsar spin geometry and illuminate its possible relation to 
SNR G343. 1—2.3. We also obtain improved estimates of the pulsar flux and nebular spectrum, 
constraining the system age and energetics. 

Subject headings: gamma rays: observations, stars: pulsars: individual B1706— 44 

1. Introduction 

PSR B1706— 44, discovered by Johnston et al. (1992) is among the most interesting pulsars for study at 
high energies. It is one of a handful of pulsars detected by EGRET in GeV 7-rays. It is quite similar to the 
Vela pulsar with a characteristic age r c = P/(2P) = 1.7 x 10 4 yr and a spindown luminosity of E 1=3 4 x 10 36 
erg/s, but is ~ 10 x more distant at d = 3e?3 kpc. Early Chandra HRC/ACIS data provided a first detection 
of X-ray pulsations and showed a compact ~ 10" surrounding pulsar wind nebula (PWN) (Gotthelf, Halpern, 
& Dodson 2002; Dodson & Golap 2002). More recent XMM-Newton spectroscopy (McGowan et al. 2004) 
has provided improved measurements of the X-ray spectrum and pulsations. Early claims that the PWN is 
detected in TeV 7-rays (Kifune et al. 1995; Chadwick et al. 1997) have not been supported by recent HESS 
observations (Aharonian et al. 2005). 

PSR B1706— 44 is superposed on a radio-bright spur of the supernova remnant G343. 1—2.3, which has a 
similar, albeit unreliable, £ — D distance of ~ 3 kpc (McAdam, Osborne & Parkinson 1993). The pulsar DM 
gives a distance of 2.3 ±0.3 kpc in the Cordes & Lazio (2002) model. Dodson & Golap (2002) have argued for 
an association. In particular, they found a faint southern extension of the SNR, which would place the pulsar 
within the full SNR boundary. They also noted an approximately N-S elongation of the X-ray PWN, pointing 
roughly back to the SNR center and argued that this would represent a trailed nebula. The required velocity 
for travel from the approximate geometric center of the SNR, about 12' away, was ~ 1000d3/r4 km/s where T4 
is the pulsar age in units of 10 4 yr. There are, however, some challenges to this SNR association. Koribalski et 
al. (1995) in an HI absorption study of the pulsar found velocity components setting lower and upper bounds 
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for the distance of d m i n = 2.4 ± 0.6 kpc and d max = 3.2 ± 0.4 kpc. However a prominent HI emission feature 
seen in the bright limb of G343. 1—2.3 at — 32km/s is not seen by Koribalski et al. in the absorption spectrum 
of the pulsar, suggesting that it lies in front of the SNR. Also scintillation studies (Nicastro, Johnston & 
Koribalski 1996; Johnston, Nicastro & Koribalski 1998) suggest a low transverse velocity for the pulsar, 
v < 89km/s. This estimate has been supported by more recent scintillation measurements (Johnston, priv. 
comm) . Thus the distances of the pulsar and SNR are still fairly uncertain. We adopt here a generic distance 
of 3 kpc in the discussion that follows, but carry through the scaling to show the distance dependence. 

Ng & Romani (2004) re-examined the AOl 50 ks HRC and ~ 15 ks ACIS-S exposures, and found that 
the compact PWN could be well modeled as an equatorial torus + polar jets, rather similar to the structures 
seen around the young Crab and Vela pulsars. They found that the brightest arc of PWN emission lay 
behind the pulsar's inferred motion from the SNR center, while a polar jet extended well in front of the 
pulsar position. This is difficult to reconcile with a bow shock/trail interpretation. On the other hand, the 
PWN symmetry axis did indeed point back to the SNR center, suggesting that a more careful evaluation of 
the connection was in order. This is particularly interesting, since comparison of the PWN symmetry (= 
pulsar spin) and proper motion axes can constrain the origin of pulsar birth kicks (Spruit & Phinney 1998; 
Lai, Chcrnoff & Cordes 2001; Romani 2004). 

We have obtained a deeper 100 ks ACIS-I exposure of PSR B1706— 44 and surroundings. The X-ray 
exposure coverage is compared to the overall geometry of G343. 1—2.3 in Figure 1. Together with new ATCA 
radio continuum imaging we are able to study the rich structure in this PWN and further constrain its 
connections with the SNR. 



2. Observations and Data Analysis 

PSR B1706-44 was observed with the Chandra ACIS-I array (4 ACIS-I chips along with the S3 and S4 
chips) on February 1-2, 2004 with standard imaging (3.2s TE) exposures. The CCD array was operated in 
'very faint' (VF) mode, allowing improved rejection of particle backgrounds. The total live-time was 98.8 ks 
and no episodes of strong background flaring were observed. Hence all data are included in our analysis. The 
pulsar was positioned near the standard aimpoint of the 13 chip and all observing conditions were normal. 
We have also compared our new exposure with the archival (February 3, 2001) 14.3 ks ACIS-S3 exposure 
(obsID 0757). As usual the backside illuminated S3 chip suffered more from particle background and after 
cutting out periods of background flares, 11 ks of clean exposure remained. All analysis was performed using 
CIAO 3.2 and CALDB 3.0.0, including automatic correction for the ACIS QE degradation. These data were 
nearly free of pile-up; the maximum pixel counts at the pulsar position indicate only 2.5% pile-up while the 
best fit model for the point source has an expected pile up fraction of ~ 3.5%. For sources with low pile-up 
we can maximize the spatial resolution of the ACIS image by removing the standard pixel randomization 
and applying an algorithm correcting the position of split pixel events (Mori ct al. 2001). This decreases 
the on- axis PSF width in our data set by > 10%. These data are compared with radio observations of the 
PWN. 
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Fig. 1. — Greyscale image of our new ACIS-I pointing of PSR BT 706-44. The contours (at 8, 10, 12 and 
14 mJy/beam) show the shell of G343. 1-2.3 from a 19 pointing 1384 MHz ATCA mosaic (Dodson & Golap 
2002). The radio map has a resolution of 70" x 47" and an rms final map noise of 0.6mJy. The X-ray PWN 
lies on a spur of radio emission. An approximate boundary of the full SNR (25' radius) and an arrow for the 
inferred PSR motion, assuming birth at the SNR center, are shown. 

2.1. Radio Imaging and Astrometry 

Data for the radio maps shown here were collected at the Australia Telescope Compact Array (ATCA) 
in Narrabri (latitude -30.3°) (Frater, Brooks & Whiteoak 1992). For the 1.4 GHz map in Figures 1 and 
2a, the data acquisition and analysis are described in Dodson & Golap (2002). For the image contours in 
Figure 3a, the data first presented in Dodson & Golap (2002) were re-imaged including the 6km baselines 
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and uniform weighting to highlight the high resolution features. The restoring beam size is 9'/0 x 7'.'8. To 
show the nebular structure, an llmJy point source PSF has been subtracted at the position of the pulsar. 
Two maxima appear flanking the pulsar position. These are unlikely to be artifacts due to pulsar variability, 
as diffractive scintillation for this pulsar is particularly weak (Johnston, Nicastro & Koribalski 1998). Since 
the data were collected in five sessions, spread over more than a year, it is in principal possible for slow 
refractive scintillation to change the pulsar flux between epochs and distort its PSF. However, each epoch 
used ~12h of integration, so any residual epoch PSF should be close to circularly symmetric, in contrast to 
the structure near the pulsar which is clearly bipolar. Further observations, with pulsar binning, have been 
requested to confirm this result. 

For the 4.8 GHz map in Figure 2b, observations were made at 4.8- and 8.6 GHz with the array in 
the standard configurations 0.75A, 1.5A and 6A on 06 Jan, 16 Feb and 11 Apr 2002. The maximum and 
minimum baselines for the 4.8-GHz data were 1 and 100 kA (angular resolutions of 3.4 to 2'.'1) for a total of 
26 hours observation. In all cases we observed the two frequencies with bandwidths of 128 MHz. We used 
the ATNF correlator mode that divides each integration's data into separate phase bins spanning the pulsar 
period. This firstly allowed the strongly pulsed point source flux to be excluded from the image and secondly 
allowed us to self-calibrate using the relatively strong point source flux from the pulsar. After data editing 
and calibrating we inverted the image with a uw-taper of 20 and deconvolved it with the full polarization 
maximum entropy task PMOSMEM in MIRIAD. 

The most important test of the SNR association would, of course, be a direct astrometric proper motion. 
With a 1.4 GHz flux of ^llmJy, PSR B1706— 44 is relatively bright. As such it is suitable for phase 
referenced VLBI astrometry, if an in-beam reference could be found. Unfortunately searches for phase 
references adequate for Australian Long Baseline Array (LBA) and US VLBA experiments have not detected 
comparison sources with compact fluxes greater than ~ lmJy. Attempts were made at external phase 
reference VLBA astrometry. However at 1.4 GHz, the nearest known reference source (2.5° away) was 
scatter broadened to ~ 50 mas. With the strong ionosphere at such low elevation, the next nearest known 
source (10° away) is too distant for effective calibration. Since the pulsar spectrum is steep, an attempt 
at VLBA astrometry at 5 GHz was also unsuccessful; at this low elevation the system temperature was 
4 — 5x nominal and only six VLBA antennae could be used, reducing the sensitivity to ~ 15% of nominal. 
So unfortunately we have only tied-array astrometry at present. Even if the pulsar does travel from the 
geometric center of G343. 1—2.3, the expected proper motion is only ~ 40mas/yr; the existing time base of 
VLA/ATCA imaging does not yet allow a serious constraint on this motion. We must conclude that a direct 
proper motion measurement awaits substantially increased (SKA or EVLA) capabilities and a long-duration, 
large base-line experiment. 



2.2. X-ray Spatial Analysis 

To show the diffuse emission surrounding PSR B1706— 44 we plot (Fig. 2a) a l-7keV image with point 
sources removed (except the pulsar) . These data are exposure corrected to minimize the chip gaps and heavily 
smoothed on a 20" scale. The diffuse emission is an edge-brightened, radius ~ 110" cavity surrounding the 
pulsar with a faint extension to the west. Contours of the 1.38 GHz radio map show good correlation with 
the radio emission in the bar crossing G343. 1—2.3 (Fig. 1). We will refer to this structure as the 'nebula'. 

Moving in toward smaller scale, in 2b we show a l-7keV image, smoothed with a 1"5 Gaussian. Point 
sources have not been removed. This shows that the cross structure fit by Ng & Romani (2004) extends 
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Fig. 2. — Left: ACIS-1 1-7 keV image with point sources (other than the pulsar) removed, exposure correction 
and 20" Gaussian kernel smoothing. Contours are from the 1.38 GHz radio map of Figure 1. Right: l'/5 
Gaussian smoothed image of the PWN with an overlay of the core of the radio nebula from a 4.8 GHz 
ATCA image (contours at 0.4, 0.8, 1.0, 1.1,.. .1.6 mJy/beam); the resolution is 20" and the image RMS is 
0.2mJy/beam. 

across ~ 1'. Narrow X-ray jets, which we refer to here as the 'outer jet' (extending south) and 'outer counter 
jet' (extending north) start ~ 10" from the pulsar and continue to ~ 30". Bracketing these is faint diffuse 
X-ray emission which we will call the 'equatorial PWN'. For comparison we draw contours of a 4.8 GHz 
ATCA image with 21" x 18" restoring beam. These observations have the pulsar 'gated out' and show that 
the radio PWN has a hollow center bracketing the equatorial PWN. Diffuse radio peaks are, in fact, seen 
just east and west of this X-ray structure. 

Finally, we show in Figure 3 a lightly smoothed image of the central region of the PWN, stretched 
to bring out the faint outer jets. The contours are drawn from a 1.38 GHz ATCA image, where the 6-km 
baselines have been weighted to produce a 9"0 x 7"8 restoring beam. A point source PSF has been subtracted 
at the pulsar position. Two local radio maxima with peak fluxes ~ 2mJy and ~ 2.5mJy bracket the 'torus' 
structure. The radio then shows a sub-luminous zone surrounding the 'equatorial PWN'; beyond ~ 30" the 
radio brightens again, as in figure 2b. No emission appears along the 'outer jets'. Indeed there appear to be 
evacuated channels in the radio emission, but improved S /N and resolution are needed to probe this sub-m Jy 
structure. The frames on the right show the innermost region of the PWN with the best-fit torus + inner 
jet model (§2.3). 

The overall geometry of the PWN is strongly reminiscent of that surrounding the Vela pulsar. In 
particular Pavlov et al. (2003) have described a series of ACIS images of the Vela nebula which show a 
torus-like structure, an inner jet and counter jet and a faint narrow outer jet system. This imaging sequence 
showed that the Vela outer jet, which is patchy and strongly bent, varies dramatically on timescales of days 
to weeks. Apparent motion of blobs within the jets suggests mildly relativistic bulk velocities and strong 
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instabilities. For PSR B1706— 44 our single sensitive image does not let us comment on variability. However 
we argue that the relatively straight and narrow jets, ~ 3x longer than those of Vela, and symmetric PWN 
structure are a consequence of a static uniform external environment and a low pulsar velocity. At 1.4- 
8.5 GHz Dodson et al. (2003) have found that the Vela PWN has two bright patches bracketing the X-ray 
torus and jets, in a structure quite similar to that in Figure 2b. Polarization imaging of the Vela radio 
structure suggests that these two patches represent the limbs of a toroidal B field structure. This implies 
that the rotation axis controls the PWN symmetry to large radii. 



Following Ng & Romani (2004) we have fitted our new ACIS image to a point source PSF, Dopplcr 
boosted equatorial torus, polar jets and uniform background. The fitting minimizes residuals using a Poisson- 
based likelihood function. Monte Carlo simulations of Poisson realizations of the best-fit model are in turn 
re-fitted to generate statistical errors and their co-variance matrices. Tabic 1 contains the best-fit values. 
The torus radius and axis inclination and position angles are r, £ and 'J', respectively. See Ng & Romani 
(2004) for the definition of the other parameters and the details of the fitting technique. In Table 1 the inner 
jet/counter-jet are constrained to lie along the torus axis in the fits. 

In addition to the statistical errors, there are certainly systematic errors, in particular induced by 



Fig. 3. — Left: Grey scale image (1'.'5 Gaussian smoothing) of the PWN, stretched to show the outer jets 
and the equatorial PWN. The contours (0.3, 0.4, 0.5, 0.6, 0.9, 1.2, 1.5mJy/beam) are from the 1.38 GHz 
ATCA image, which has a map rms of 0.2 mJy/beam and a resolution of 9'.'0 x 7'.'8. The pulsar point source 
has been subtracted. Two peaks of radio emission bracket the torus. These lie in a cavity which surrounds 
the equatorial PWN and jets. Middle: grey scale image of the inner PWN (above) and of the best-fit point 
sourcc+torus+ (inner )jets model (Table 1) to the same scale(below). At the far right, the residual image is 
shown, with contours indicating the excess counts above (North of) the torus. 



2.3. 



Nebula Structure Fits 
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Tablc 1: Torus Fit Parameters with la Statistical Errors 

* (deg) C (deg) r (arcsec) 5 /3 Point Source 1 ^ Torusf Jetf Counter-jetf 

163.6 ± 0.7 53.3±i;| 3.3^'^ L0* 0.70 ±0.01 2897+fl 1221+^ 185+^ 325±fl 
* - held fixed f counts 



unmodclcd PWN components. For example, it is clear that there are counts in excess of the torus+jet model 
in a cap surrounding the inner counter-jet. Interestingly similar structure is seen in the Crab PWN. We have 
made an attempt to constrain the systematic biases by modifying the fitting model. For example, allowing 
the (inner) jet and counter-jet to have a free position and amplitude shifts the best-fit position angle to 
= 165 ± 0.5° and the inclination to £ = 56.7 ± 1.0°. We therefore infer systematic errors about 3x larger 
than our rather small statistical errors. 

We have also measured the outer jet /counter- jet system. Minimizing the residual to a 1-D line passing 
through the pulsar, the two jets together lie at Pouter = 169.4 ± 0.15°. If the jets are fitted separately, we 
obtain * outer = 168.4 ± 0.2° and 170.9 ± 0.2° for the outer jet and counter-jet, respectively. Thus the two 
jets are mis-aligned at the ~ 8a level. A fit to the count distribution about the best fit axis shows that the 
narrow outer counter-jet has a Gaussian FWHM across the jet of 2'.'3 ± 0'.'2. The outer jet appears broader 
at the base with an initial width of 4'/9 ± 0'.'5, continuing at FWHM=2"7 ± 0"3 for its outer half. These 
estimates have been corrected for the telescope PSF, which is quite uniform this close to the aimpoint. 

It is important to note that at the observation roll angle, the read-out direction lies at ^ = 168.7°! Due 
consideration, however, shows that the jet structure cannot be produced by the read-out trail. First the jets 
cover only ~ 1'; the read-out excess should cover the full 13 chip. Second the pulsar provides only ~ 2900 
1 — 7keV counts. The read-out trail (out-of-time) image of this source should contribute only 36 counts over 
the full 8.3' strip across 13 and ~ 2.5 counts in the 'jet regions'; the outer jet and counter-jet have 92 and 
93 1 — 7kcV counts, respectively. Finally the outer jets are much harder than the soft X-ray emission from 
the pulsar; indeed with a mean detected photon energy ~ 2 keV these are the hardest extended features in 
the image. 

The over-all system, showing an asymmetric torus, broad inner jets and narrow outer jets is, of course, 
very similar to the Vela PWN as studied with Chandra by Pavlov et al. (2003). We will discuss the comparison 
with the Vela system in §3, highlighting the differences. We interpret these as suggesting that the PSR 
B1706— 44 PWN has developed from a low velocity pulsar. 



2.4. Spectral Analysis 

For the best possible constraints on the source spectrum, we have reprocessed both the ~ llks cleaned 
ACIS-S data set and our new ~ 100 ks ACIS-I data set with the new time-dependent gain adjustment and 
CTI correction available in CI AO 3.2. The updated RMFs should in particular improve the low energy 
calibration, important for obtaining the best estimates of Nh- As noted above, in these data sets the pile-up 
was negligible at ~ 3%. To model the aperture corrections, 10 PSFs with monochromatic energies from 0.5 
to 9.5keV were simulated using the Chandra Ray- Tracer program, ChaRT. The enclosed energy fraction as 
a function of radius was fitted to a linear function of energy and this was used to correct the ARFs used in 
the spectral fit. In extracting the pulsar spectrum, an aperture of radius 1" was used to minimize nebular 
contamination. Results from the combined fits of the ACIS-I and ACIS-S pulsar data sets are listed in 
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Table 1; the spectral fits are substantially better for composite models with both thermal and power-law 
components. Spectral parameter errors are projected multi-dimensional la values. We quote both absorbed 
and unabsorbed fluxes. As is often the case with low statistics X-ray spectra, projected (multi-dimensional) 
errors on the fluxes are very large due to spectral parameter uncertainties. Thus, we follow other authors in 
quoting flux errors as la single parameter values. 

Table 2: Spectral Fits to PSR Bf706-44 



Model 


N„ 


r 


abs.flux 


unabs. flux 


rpOO 


R°° 


abs. flux 


unabs. flux 


X 2 /dof 




10 21 cm~ 2 




JO. 5-8 


JO. 5-8 


MK 


km 


Jo. 5-8 


Jo. 5-8 




PL 


3.8 ± 0.3 


3.061!!;- 


2.01 ± 0.04 


4.74 ± 0.1 










145.9/99 


BB 


1.0 ± 0.3 








3.7 ± 0.2 


u '*-0.11 


1.42 ± 0.03 


1.76 ± 0.04 


288.4/99 


PL+BB 


4 y+0.8 
*■ ' -0.7 


1.62 ± 0.2 


1.70 ± 0.07 


2.25 ± 0.1 




2 4+ 1 ' 6 


0.83 ± 0.03 


3.85 ± 0.12 


40.6/97 




5.0* 


1.68 ± 0.17 


1.70 ± 0.07 


2.33 ± 0.1 


1.93 ± 0.06 


o O + 0.76 
Z ' 5 -0.70 


0.82 ± 0.03 


4.40 ± 0.14 


40.8/98 


PL+Atm 


5 9+ 09 

°' y -0.8 


1.67 ± 0.2 


1.67 ± 0.07 


2.34 ± 0.1 


7Q + 013 


13.1* 


0.85 ± 0.03 


6.8 ± 0.2 


40.7/97 




5.0* 


1.56 ± 0.17 


1.67 ± 0.07 


2.2 ± 0.1 


0.91 ± 0.05 


13.1* 


0.88 ± 0.03 


4.92 ± 0.15 


40.8/98 



t 0.5-8 keV fluxes in units of 10 13 erg/cm 2 /s * held fixed 



To get the best constraints on the point source spectrum, Table 2 gives fits with Njj held to the 
value from the power law fits to the extended emission. We have also compared our results with the XMM- 
Newton fitting of McGowan et al. (2004), by fitting counts in the 20" aperture used in that observation. 
Our parameters and fluxes for the thermal component are generally in very good agreement. However, since 
this aperture contains much of the torus and central PWN, XMM-Newton substantially overestimates the 
non-thermal flux for the point source. Their fit power law flux corresponds to 7.5 x f _13 erg/cm 2 /s (0.5- 
8keV, unabsorbed). In the 20" aperture we find 9.6 x 10 _13 erg/cm 2 /s (0.5-8keV); while the small Chandra 
point source aperture gives 2.3 x 10 _13 erg/cm 2 /s (0.5-8keV) for the power law component. We find a similar 
~ 3x excess in the power-law + atmosphere flux for the fit to the large XMM-Newton aperture. Conversion 
of the power law flux observed in our small point source aperture to the XMM-Newton band shows that 
the expected PN+MOS (0.2-10 keV) count rate is 18% of the total (power law + thermal) counts in the 20" 
aperture. In soft (0.2-1.35 keV) and hard (1.35-10kcV) bands the predicted fraction of the counts from the 
power law are 12% and 21%, respectively. However the light curves of McGowan et al. (2004) show that the 
pulse fractions are 21% (soft), 12% (hard) and 11% (total). Since the small aperture power law produces 
only 12% of the soft counts but 21% are pulsed, there must a thermal pulse component. Conversely, since 
the power law produces 21% of the hard band flux, but this only has a pulse fraction of 12%, some of the 
power law counts must be unpulsed. Extrapolation of the PWN count excess above the point source PSF 
in the sub-luminous zone at 2—3" produces < 1% of the point source aperture counts. Thus the larger scale 
torus emission does not contribute significantly to the point source power law and cannot account for its 
unpulsed component. This suggests that part of the magnctosphcric emission is nearly isotropic or that 
there is a very compact (< 1") PWN component at the pulsar position. 

For the thermal component, the fit flux gives an emitting area (effective radius) as a function of distance. 
Our fit to a pure blackbody gives R e ff — 2.8e? 3 km. Thus for reasonable distances, this flux represents hot 
T ~ 2 x 10 6 K emission from a small fraction of the stellar surface (~ 4.5% for an R°° = 13.1km star). 
The light element neutron star atmosphere models, such as the pure H 10 12 G model grid used here (Zavlin 
et al. 1996), have large Wien excesses. When fit they give lower T e ff. Also, the black body departures 
allow one, in principle, to fit both the surface redshift and radius. In practice, these are typically highly 
degenerate in CCD-quality data. We assume here a generic surface radius of R s = 10 km, corresponding 
to R x = R s (l - 2GM/R.SC 2 )- 1 / 2 = 13.1km. With N H free (giving 5.9 ± 0.9 x 10 21 c m - 2 ) our thermal 
flux normalization gives a radiating radius of Roo — 27.4d 3 km, which is difficult to reconcile with expected 
neutron star radii for any d > 1.8 kpc. However, when Nh is fixed at the nebular value of 5 x 10 21 cm~ 2 , wc 
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get an effective radius of = I6.lt/3km, which is tolerable even at our nominal 3kpc distance. 

Analysis of the low signal-to-noise, extended flux depends critically on the background subtraction. 
Given the limited statistics, only simple absorbed power-law fits were attempted for all non-thermal sources. 
The results are listed in Table 3. For consistency, all fits are to the 0.5-8 keV range and we quote both 
absorbed and unabsorbed fluxes. 



Table 3: Spectral fits to extended sources 
Object N H T 

10 21 cm- 2 



abs. flux unabs. flux x 2 /dof 



./r 



+0.00 

-0.08 



Nebula 



.OiO.4* 1.77 



Equatorial PWN 5.0 ± 0.4* 1.57 ±0.08 
Torus 5.0 ±0.4* 1.48 ±0.08 



5.6 ±0.23 

2.4 ±0.07 

1.5 ±0.05 



+0.14 



7.9 ± 0.32 67.4/77 
3.2 ± 0.09 70.3/81 
1.9 ±0.07 39.8/53 



-0.13 



Jet 5.0* 1.26 

Counter Jet 5.0* 1.39 ±0.10 

Outer Jets 5.0* 1.26 ±0.18 

* 0.5-8 keV fluxes in units of 10~ 13 erg/cm 2 /s 



0.42 ±0.03 0.52 ±0.04 12.0/18 
0.81 ±0.04 1.0 ±0.05 10.8/18 
0.27 ±0.03 0.33 ±0.03 11.9/18 

* held fixed * Simultaneous Nh fit. 



Note that there is significant softening of the extended emission as one progresses to larger scales and 
that the jet components appear to be the hardest of all. This is certainly consistent with the idea that 
the central pulsar supplies fresh energetic electrons and that synchrotron burn-off increasingly softens the 
spectrum as older populations are viewed in the outer PWN. Again this trend is common in the well- 
measured young PWNe. Allowing the photon index to vary for the different nebula components, the best 
fit to a global absorption value for the extended emission gives us our fiducial Nh — 5 x 10 21 cm~ 2 . This is 
consistent with free-fit values for the point source, but given complexities of the composite thermal±power 
law model, we consider the nebular fit value more robust. Note that with DM=75.7 cm~ 3 pc, the H/n e « 21 
for this sight-line is large, but not unprecedented for low |6| pulsars. This is also consistent with the Hi 
absorption measurements and a fiducial SNR distance ~ 3 kpc, given the appreciable uncertainties. 



3. Interpretation and Conclusions 

A number of authors have discussed the evolution of a PWN within an expanding supernova remnant. 
For example, van der Swaluw (2001) and Chevalier (2005) describe the early evolution when the supernova 
ejecta are in free expansion. Later, after the remnant interior is heated by the passage of the reverse shock, 
the PWN evolves within the Sedov phase supernova remnant whose radius is i?sNR = 1.17(Eo/ p) 1 / 5 ^/ 5 for 
an explosion energy Eq in a 7 = 5/3 medium of density p. PSR B1706— 44 has a characteristic age 10 4 T4 yr 
with r 4 w 1.7, so G343. 1—2.3 should be safely in the Sedov phase with an expected angular size 

0snr « 16'(£ 5 i/no) 1/5 *4 /5 M (1) 

for a supernova releasing energy E — 10 51 i^5ierg in an external medium density n cm~ 3 , at a true age 
10 4 £4y at a distance 3^3 kpc. The observed size then implies E51 llnot^ 2 d^, requiring a fairly energetic 
explosion for d > 2 kpc. During the Sedov phase the interior pressure is 

P SNR « lO- 9 ^ V V /5 g/cm/s 2 (2) 

and is relatively constant away from the SNR limb. 
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The pulsar blows a wind bubble within this SNR interior, whose radius is Ppwn ~ (-E'*/-E'o) 1 / 3 -Rsnr for 
a PWN bubble energy E* = fEr c (van der Swaluw & Wu 2001). Although the accuracy of this dependence 
of PWN radius on pulsar injection energy has been questioned (Blondin, Chevalier & Frierson 2001), wc 
adopt it for the following estimates. With the observed ratio of radii, Ppwn/Psnr = 1.8'/25', we obtain 
P* = 3.7 x 10~ 4 Po, i.e. this PWN has quite low internal energy. This is also reflected in the low radio and 
X-ray fluxes. Together we use these estimates, the observed size of the SNR, equation (1) and the measured 
P 36 = 3.4 and r c = 1.75 x 10 4 y to write / « (Rpwn/Rsnr) 3 E /(Et c ) = lln^dl,. Now, if the PWN is 
adiabatic and we assume spindown with constant B and braking index n — 3 from an initial period P , wc 
find that the total energy in the plerion is [(P/Po) 2 — 1]Et c . Then, setting / = (P/Pq) 2 — 1 and eliminating 
the true age t using t = r c [l — (Po/P) 2 ] for magnetic dipolc spindown, wc obtain a constraint on the initial 
spin period 

[1 - (P /P) 2 ] 3 /(P)/P) 2 = 0.68 n 4 (3) 

which has a solution of P = 0.61P = 62 ms for d=3kpc, and P = 0.79P = 80 ms for d=2kpc. The 
corresponding true ages are 0.61t c (1.1 x 10 4 y) and 0.38r c (0.67 x 10 4 y), respectively. These numerical 
values are for no = 1 and the density dependence from Equation (3) is quite weak, van der Swaluw & Wu 
(2001) present a similar sum for P , assuming a known E ; the above formulation emphasizes the sensitivity 
to the poorly known d. Note that with the large implied initial period, the integrated PWN energy is quite 
comparable to the present spin energy, with / w 1.7 at d=3kpc and / w 0.62 at d=2kpc. So the spindown 
luminosity is roughly constant in the adiabatic phase and the PWN growth is closer to t 11 / 15 than to the 
i 3 / 10 law appropriate for impulsive energy injection (van der Swaluw 2001). 

Inside this wind bubble, the Sedov interior pressure confines the PWN, giving rise to a termination 
shock at 

9 WS « (P/47rcP SNR ) 1 / 2 /d. (4) 

which results in 6*ws ~ ^''^E^g 2 E 5 ^ 5 n 3 ^ 1 °t^ 5 d 3 ~ 1 . If we apply the SNR estimate for Eq above, this becomes 
6> ws « 0'. / 72El( i 2 ng 1/2 t i d^ 2 . Then, using P36 « 4 and applying the age estimate following Equation (3) we 
get 6*ws ~ l'.'5n ^ 2 (d=3kpc) or 6*ws ~ 2'/ln ^ 2 (d=2kpc). These estimates are reasonably consistent 
with the observed 3" torus radius, especially since an cquatorially concentrated flow should have a stand-off 
distance 1.5-2 x this spherical scale. The polar jets can have an initial shock at somewhat larger angle, with 
the resulting pitch angle scattering illuminating the jets somewhat further from the pulsar. 

Of course, this bubble is offset from the center of G343. 1—2.3 at R — 0.5 Rsnr (figure 1). This is inside 
the ~ 0.68Psat.r where van der Swaluw, Downes & Keegan (2004) note that the increasing density causes 
the pulsar to be supersonic, so a bow shock should not have yet formed. These authors however compute 
numerical models of a fast moving pulsar in a SNR interior. As the pulsar moves, the PWN should become 
highly asymmetric with a 'relic PWN' at the SNR center and the pulsar placed near the leading edge of the 
PWN; see van der Swaluw, Downes & Keegan (2004) figures 7 and 8. We see no PWN structure near the 
geometric center of G343. 1—2.3 and, if the 'bubble nebula' is identified with the shocked pulsar wind, the 
pulsar is certainly not offset from its center along the proper motion axis (away from the SNR center). So 
these models are an inadequate description of G343. 1—2.3. From Figure 2, the pulsar is well centered in the 
bubble nebula, with any offset from its center along the axis to the SNR substantially less than 30A 30 arcsec. 
Thus 

v < 40A 30 d 3 /U km/s, (5) 

and the pulsar cannot have moved far from the explosion center. This is, of course, consistent with the 
scintillation results. 
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We can reconcile the symmetric PWN with the offset SNR shell if we assume that the pulsar progenitor 
exploded toward the edge of a quasi-spherical cavity. One scenario (also posited by Gvaramadze 2002, 
Bock & Gvaramadze 2002) that can associate the low velocity pulsar with G343. 1—2.3 is to assume that 
the progenitor star had a stellar wind of mass loss rate M_8lCr 8 M Q /yr and wind speed 10 8 w w s cm/s over 
t ~ til0 7 y, typical of the <~ 10M Q stars that dominate the pulsar progenitors (Maeder 1981). This evacuates 
a stellar wind bubble of size 

SW = 46' (M_ 8 ^ 8 /n ) * " t 3 7 /5 /d 3 . (6) 

During the main sequence lifetime, the star moving at lOwgkm/s travels ~ 2°VQtT/d 3 and so it can easily 
traverse its wind bubble. Thus, one can imagine an off-center supernova in a nearly symmetric stellar wind 
bubble of radius ~ 25': the supernova blast wave expands to fill the bubble, passing to the Sedov phase 
near its present radius. The supernova produces a neutron star with little or no kick, placing the pulsar 
near its present position. This has the added advantage of accommodating the rather large SNR size with a 
more modest energy of a fewxlO 51 erg. The PWN energy and size estimates above would then be somewhat 
amended; this would require a careful numerical simulation. 

For the reasons detailed in the introduction, it is not yet clear that PSR B1706— 44 and G343. 1—2.3 are 
associated. So for completeness we can consider the case when the shocked pulsar wind blows an adiabatic 
bubble in a static, low P oxt external medium (Castor, McCray & Weaver 1975). If we assume that the pulsar 
was born (sans SNR) or entered a confining region of the ISM ~ 10 4 y ago and that since then it has been 
spinning down at the present energy loss rate, we find that it will blow a bubble of angular size 

0bn » 0.76(Et 3 / P ) 1/5 /d « 120"(E 36 /n Q )^ 5 T 3/5 /d 3 . (7) 

These estimates change somewhat for a pulsar born at Pq <C P; since we are not making the association 
with the SNR G343. 1—2.3, we can make no estimate of the initial spin period. As first noted by Dodson 
& Golap (2002), the ~ 4' wide radio spur across the face of G343. 1—2.3 has the approximate scale of such 
a 'bubble nebula'. If the PWN stays unmixed (relativistic) then the interior of the bubble will have a 

/ \l/5 

pressure Pbn ~ Et/(4irR 3 ) w 1.6 x 10~ 10 (nQ-Efg/rf ) g/cm/s 2 . In turn, the torus termination shock in 
this medium is at 

9 WS « 2'!9(E 36 /n )V W T% /5 /d 3 . (8) 

If (e.g. through Rayleigh- Taylor instabilities) the pulsar wind is well mixed with the swept up gas, the 
adiabatic thermal pressure would be ~ 2x larger. Interestingly, the angular scales for this scenario are also 
reasonably compatible with the observed torus and bubble nebula size. Of course this scenario leaves open 
the question of the pulsar origin. Again the pulsar would need to have a quite low velocity to produced the 
observed symmetry. 

Turning to the spectral results, we note that Possenti et al. (2002) fit a correlation between spindown 
energy and the PSR+PWN luminosity: L(2 - lOkcV) = 1.8 x 10 38 £4 34 erg/s. For the PSR B1706-44 
parameters this predicts a flux _ f (2 — 10 keV) = 4.7x 10~ 12 <i3~ 2 erg/cm 2 /s. The observed 2-10 keV flux is in fact 
~ 1.7 x 10~ 12 erg/cm 2 /s, even including the outer 'bubble nebula'; without this component it is half as large. 
These correlations are not very accurate, but this does imply that the PSR B1706— 44 PWN is substantially 
under-luminous for any distance less than 3 kpc. Gotthelf (2003) has derived correlations between the 
pulsar spindown power and the pulsar/PWN spectral indices. His relation predicts TpgR, = 0.63 ± 0.17 
(substantially smaller than our power law index ~ 1.6 ± 0.2) and Tp-yvN = 1-3 ± 0.3 (not inconsistent with 
the values measured for the torus and equatorial PWN). 
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As described in §2.4, the spectrum softens appreciably from the central torus to the outer bubble nebula. 
This suggests increased aging of the synchrotron population. Figures 2 and 3 show that the bulk of the radio 
emission lies in the 'bubble nebula' region. So we can take the radio flux and spectral index from Giacani et 
al. (2001) and compare with our nebula X-ray flux (Fig. 4). Comparing the radio spectral index an — 0.3 
with the best fit X-ray index ax — 0.77, shows a break quite close to the Aa = 0.5 expected from synchrotron 
cooling. The extrapolated intersection of these power laws gives a break frequency of Log[^s(Hz)] = 12.2+^i- 
For the fiducial pulsar age of ~ 1.7 x 10 4 y, this corresponds to a nebula field of 1.4+ ; 6 x 10~ 4 G. Note that 
the magnetic pressure from this (photon flux-weighted) average field is ~ 8 x 10 _10 g/cm/s 2 , somewhat 
larger than the nebula pressure estimated from its radius. This may indicate field compression in the nebula 
limb. In general, if the mean nebula field is 10~ 4 _B_4 G for a nebula of angular radius <~ 100"6>i O o, the 
total nebula field energy is Eb ~ 1.5 x 10 47 _B 2 4 (6*iooe?3) 3 erg. This is comfortably less than the present spin 
energy £psR ~ 2 x 10 48 erg, so the nebula can be easily powered even if the pulsar was born close to its 
present spin period. We find that this cooling break field is substantially larger than the equipartition field 
of 10 — 15/iG inferred for the radio and X-ray emitting populations (also the minimum equipartition nebula 
energy ~ 9 x 10 45 erg is substantially smaller). The cooling break field can also be compared to that expected 
from simple radial evolution of the pulsar surface field: if this field B* = 3 x 10 12 G falls off as r -3 to the 
light cylinder, then as 1/r to the wind shock where it is compressed we get -Bws ~ 3S*r 3 /(r 2 jC rws) ~ lmG. 
If it continues to fall off as 1/r beyond this we get a field at the limb of the bubble nebula of ~ 30^G. So 
the best we can do is to infer a mean nebular field ~ 10 — 30x the equipartition value, with some generation 
of new field beyond the torus wind shock. The energetic requirements for this field, required to match the 
vb cooling break, are comfortably less than the energy available from PSR B1706— 44. These field estimates 
are consistent with the non-detection of TeV ICS flux from this source (Aharonian et al. 2005). 



The narrow outer jets also have a power-law spectrum and are almost certainly synchrotron-emitting. 
For a reasonable 0.1i3_4mG field, the observed X-rays of E 1 = 1.5_B_4r 2 5 kcV require substantial 
energies, with T e = 3 x 10 7 r7.5 near the radiation-reaction limited primary Lorentz factor inferred for 
many polar cap (Muslimov & Harding 2003) and outer magnctosphere (Romani 1996) pulsar models. Since 
the jet is narrow, confinement of these pairs imposes a (not very restrictive) lower bound on the jet field 

1/3 

B_ 4 > 0.075E^°/(d 3 6 w ) where the maximum observed jet photon has E ~ 5E$ keV and the observed outer 
jet half-width is ^^arcsec. A more restrictive upper limit on the mean jet field comes from the observation 
that the jets do not soften noticeably before their end ~ 30$3oarcsec from the pulsar. If we assume a jet 
bulk speed /3c, then arguing that the flow time is shorter than the synchrotron cooling time gives us the 

2 / 3 —1/3 



limit B_ 4 < 8.5 (J^) E 5 



When the observed jet spectrum has an energy index of a = T — 1 « 0.3, we infer a power-law spectrum 
of in the jet N(T e )dT e = KT~ p dT e , withp = 2a+l « 1.6. We can then make an estimate of the minimum 
jet luminosity, i.e. at 'equipartition' when B 2 = 6irm e c 2 J T e N(T e )dT e . Given the observed combined outer 
jet luminosity (0.5-8keV) L = 47rrf 2 / oj w 2.7 x 10 31 erg/s and emitting volume V w 2 x 8 L x ir8 w d 3 « 
1.1 x 10 52 (#t,/20)# 2 d 3 cm 3 with the angles in arcsec, we can estimate the equipartition field for an isotropic 
plasma as 

"l87r (2irm e c\ 1 l 2 2(1 - a) 1 - (E min / E^ 1 - 2 ^' 2 L 1 ' 



B eq 



(9) 



IT V ^max J 1 -2a 1 - (£ m i„/£ m ax) (1 - a) V 

where the observed photon spectrum runs from E min to E max . For the observed flux this gives 

B cq w 0.25 x 10- 4 <?(a)[^ 2 d3/20]- 2 / 7 G, (10) 



- 13 - 



where q(a = 0.3) — 1 is a weak function of a. The corresponding minimum energy flux for the outer jet is 

L oj - 8 x lO 33 /34 2/7 ^°/ 7 (0 L /2O)- 2 / 7 erg/s, (11) 

where the jet bulk velocity is /3c. This is ~ 10~ 3 E per jet and will, of course, be larger if the jet flow includes 
ions. Interestingly, if the pulsar couples roughly isotropically to the PWN, then the corresponding fraction 
of the outflow should subtend a half angle of ~ 5°. This is somewhat smaller than the angle subtended by 
the inner jets, but ~ 3x larger than the ~ 1" width of the ends of the jet - there is substantial collimation 
of the jet energy flux. 

We have argued that a low PSR velocity can explain the symmetry of the PWN. The central location 
of the pulsar and spherical post-shock flow may also allow the equatorial toroidal structure and polar jets 
to propagate undisturbed to large radii. We do, however, measure a small misalignment of the outer jets, 
corresponding to a deflection of 6>d c = 1.3 ± 0.15° for each. If we imagine a pressure acting along the jet's 
~ 30" length, then the required perturbation is 5P w L jtan6*dc/ (/3cA>j) w 5 x 10^ 14 L 34 / (f39 30 w dl)g/cm/s 2 , 
where A j is the jet's cross sectional area. This is only ~ 10~ 3 of the total pressure in the nebula. It could 

]/2 

be due to ram pressure if the shocked nebular medium flows to the west at v <~ 1.7n neb ' km/s. 

Our X-ray measurements have established the PWN symmetry axis, presumably reflecting the pulsar 
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Fig. 4. — Outer PWN ('bubble nebula') spectral energy distribution (SED). Radio data are from Giacani 
et al. (2001), the TeV upper limit is from Aharonian et al. (2005). The spectra of the inner, younger PWN 
components are plotted for comparison. The extension of the radio PL and the best fit nebula (outer PWN) 
PL meet at vb ~ 1.6 x 10 12 Hz; the indices are consistent with a Aa = 0.5 cooling break. 
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spin axis, to very high precision. Unfortunately our original goal of relating this to the proper motion axis 
remains unfulfilled. It is true that the torus symmetry axis points roughly toward the center of G343. 1—2.3, 
confirming the estimates from earlier Chandra data. However, the PWN symmetry about the pulsar and the 
low scintillation velocity suggest a very low transverse speed < 40km/s, which would preclude a birth site as 
distant as the SNR center. This low speed makes a direct proper motion challenging, but allows latitudinal 
asymmetries in the PWN flow to propagate undisturbed to fairly large radius, where they can be imaged 
with Chandra. Thus, study of this PWN offers some good opportunities to probe outflow dynamics and jet 
collimation. Study of this, and similar, PWNe may prove useful electrodynamic analogs of the 10 6 x more 
powerful AGN jets. A viable scenario for maintaining the G343.1— 2.3/PSR B1706— 44 association posits 
a supernova event near the present pulsar site, with the remnant inflating a pre-existing off-center cavity. 
However, the residual (small) proper motion could then have any direction. In fact, the faint extension of 
the PWN (bubble nebula) to the west and the increased radio surface brightness to the east might suggest a 
rather slow pulsar motion at PA ~ 80°. This would be nearly orthogonal to the torus symmetry axis. With 
the large Pq estimated here, this could be construed as suggesting poor rotational averaging of a birth kick 
(Ng & Romani 2004). So the PWN/SNR geometry offers both aligned and orthogonal axes. Only a sensitive 
astrometric campaign can detect or limit the pulsar motion and resolve this ambiguity. 

This work was supported in part by NASA grants SAO G04-5060X and NAG5-13344. We thank the 
referee for several careful readings, which resulted in substantial changes. We also thank R.N. Manchester 
who provided a current pulsar ephemeris for the VLBI experiment. 
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